Identification of direct residue contacts in protein-protein interaction by message passing

Understanding the molecular determinants of specificity in protein-protein interaction is an outstanding challenge of postgenome biology. The availability of large protein databases generated from sequences of hundreds of bacterial genomes enables various statistical approaches to this problem. In this context covariance-based methods have been used to identify correlation between amino acid positions in interacting proteins. However, these methods have an important shortcoming, in that they cannot distinguish between directly and indirectly correlated residues. We developed a method that combines covariance analysis with global inference analysis, adopted from use in statistical physics. Applied to a set of >2,500 representatives of the bacterial two-component signal transduction system, the combination of covariance with global inference successfully and robustly identified residue pairs that are proximal in space without resorting to ad hoc tuning parameters, both for heterointeractions between sensor kinase (SK) and response regulator (RR) proteins and for homointeractions between RR proteins. The spectacular success of this approach illustrates the effectiveness of the global inference approach in identifying direct interaction based on sequence information alone. We expect this method to be applicable soon to interaction surfaces between proteins present in only 1 copy per genome as the number of sequenced genomes continues to expand. Use of this method could significantly increase the potential targets for therapeutic intervention, shed light on the mechanism of protein-protein interaction, and establish the foundation for the accurate prediction of interacting protein partners.Comment: Supplementary information available on http://www.pnas.org/content/106/1/67.abstrac

Purification and Activity Testing of the Full-Length YycFGHI Proteins of Staphylococcus aureus

Background: The YycFG two-component regulatory system (TCS) of Staphylococcus aureus represents the only essential TCS that is almost ubiquitously distributed in Gram-positive bacteria with a low G+C-content. YycG (WalK/VicK) is a sensor histidine-kinase and YycF (WalR/VicR) is the cognate response regulator. Both proteins play an important role in the biosynthesis of the cell envelope and mutations in these proteins have been involved in development of vancomycin and daptomycin resistance. Methodology/Principal Findings: Here we present high yield expression and purification of the full-length YycG and YycF proteins as well as of the auxiliary proteins YycH and YycI of Staphylococcus aureus. Activity tests of the YycG kinase and a mutated version, that harbours an Y306N exchange in its cytoplasmic PAS domain, in a detergent-micelle-model and a phosholipid-liposome-model showed kinase activity (autophosphorylation and phosphoryl group transfer to YycF) only in the presence of elevated concentrations of alkali salts. A direct comparison of the activity of the kinases in the liposomemodel indicated a higher activity of the mutated YycG kinase. Further experiments indicated that YycG responds to fluidity changes in its microenvironment. Conclusions/Significance: The combination of high yield expression, purification and activity testing of membrane and membrane-associated proteins provides an excellent experimental basis for further protein-protein interaction studies an

Feedback control architecture and the bacterial chemotaxis network.

PMCID: PMC3088647This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Bacteria move towards favourable and away from toxic environments by changing their swimming pattern. This response is regulated by the chemotaxis signalling pathway, which has an important feature: it uses feedback to 'reset' (adapt) the bacterial sensing ability, which allows the bacteria to sense a range of background environmental changes. The role of this feedback has been studied extensively in the simple chemotaxis pathway of Escherichia coli. However it has been recently found that the majority of bacteria have multiple chemotaxis homologues of the E. coli proteins, resulting in more complex pathways. In this paper we investigate the configuration and role of feedback in Rhodobacter sphaeroides, a bacterium containing multiple homologues of the chemotaxis proteins found in E. coli. Multiple proteins could produce different possible feedback configurations, each having different chemotactic performance qualities and levels of robustness to variations and uncertainties in biological parameters and to intracellular noise. We develop four models corresponding to different feedback configurations. Using a series of carefully designed experiments we discriminate between these models and invalidate three of them. When these models are examined in terms of robustness to noise and parametric uncertainties, we find that the non-invalidated model is superior to the others. Moreover, it has a 'cascade control' feedback architecture which is used extensively in engineering to improve system performance, including robustness. Given that the majority of bacteria are known to have multiple chemotaxis pathways, in this paper we show that some feedback architectures allow them to have better performance than others. In particular, cascade control may be an important feature in achieving robust functionality in more complex signalling pathways and in improving their performance

Coevolved mutations reveal distinct architectures for two core proteins in the bacterial flagellar motor

Switching of bacterial flagellar rotation is caused by large domain movements of the FliG protein triggered by binding of the signal protein CheY to FliM. FliG and FliM form adjacent multi-subunit arrays within the basal body C-ring. The movements alter the interaction of the FliG C-terminal (FliGC) "torque" helix with the stator complexes. Atomic models based on the Salmonella entrovar C-ring electron microscopy reconstruction have implications for switching, but lack consensus on the relative locations of the FliG armadillo (ARM) domains (amino-terminal (FliGN), middle (FliGM) and FliGC) as well as changes during chemotaxis. The generality of the Salmonella model is challenged by the variation in motor morphology and response between species. We studied coevolved residue mutations to determine the unifying elements of switch architecture. Residue interactions, measured by their coevolution, were formalized as a network, guided by structural data. Our measurements reveal a common design with dedicated switch and motor modules. The FliM middle domain (FliMM) has extensive connectivity most simply explained by conserved intra and inter-subunit contacts. In contrast, FliG has patchy, complex architecture. Conserved structural motifs form interacting nodes in the coevolution network that wire FliMM to the FliGC C-terminal, four-helix motor module (C3-6). FliG C3-6 coevolution is organized around the torque helix, differently from other ARM domains. The nodes form separated, surface-proximal patches that are targeted by deleterious mutations as in other allosteric systems. The dominant node is formed by the EHPQ motif at the FliMMFliGM contact interface and adjacent helix residues at a central location within FliGM. The node interacts with nodes in the N-terminal FliGc α-helix triad (ARM-C) and FliGN. ARM-C, separated from C3-6 by the MFVF motif, has poor intra-network connectivity consistent with its variable orientation revealed by structural data. ARM-C could be the convertor element that provides mechanistic and species diversity.JK was supported by Medical Research Council grant U117581331. SK was supported by seed funds from Lahore University of Managment Sciences (LUMS) and the Molecular Biology Consortium

Evolution of Taxis Responses in Virtual Bacteria: Non-Adaptive Dynamics

Bacteria are able to sense and respond to a variety of external stimuli, with responses that vary from stimuli to stimuli and from species to species. The best-understood is chemotaxis in the model organism Escherichia coli, where the dynamics and the structure of the underlying pathway are well characterised. It is not clear, however, how well this detailed knowledge applies to mechanisms mediating responses to other stimuli or to pathways in other species. Furthermore, there is increasing experimental evidence that bacteria integrate responses from different stimuli to generate a coherent taxis response. We currently lack a full understanding of the different pathway structures and dynamics and how this integration is achieved. In order to explore different pathway structures and dynamics that can underlie taxis responses in bacteria, we perform a computational simulation of the evolution of taxis. This approach starts with a population of virtual bacteria that move in a virtual environment based on the dynamics of the simple biochemical pathways they harbour. As mutations lead to changes in pathway structure and dynamics, bacteria better able to localise with favourable conditions gain a selective advantage. We find that a certain dynamics evolves consistently under different model assumptions and environments. These dynamics, which we call non-adaptive dynamics, directly couple tumbling probability of the cell to increasing stimuli. Dynamics that are adaptive under a wide range of conditions, as seen in the chemotaxis pathway of E. coli, do not evolve in these evolutionary simulations. However, we find that stimulus scarcity and fluctuations during evolution results in complex pathway dynamics that result both in adaptive and non-adaptive dynamics depending on basal stimuli levels. Further analyses of evolved pathway structures show that effective taxis dynamics can be mediated with as few as two components. The non-adaptive dynamics mediating taxis responses provide an explanation for experimental observations made in mutant strains of E. coli and in wild-type Rhodobacter sphaeroides that could not be explained with standard models. We speculate that such dynamics exist in other bacteria as well and play a role linking the metabolic state of the cell and the taxis response. The simplicity of mechanisms mediating such dynamics makes them a candidate precursor of more complex taxis responses involving adaptation. This study suggests a strong link between stimulus conditions during evolution and evolved pathway dynamics. When evolution was simulated under conditions of scarce and fluctuating stimulus conditions, the evolved pathway contained features of both adaptive and non-adaptive dynamics, suggesting that these two types of dynamics can have different advantages under distinct environmental circumstances

Genetic Analysis of the Functions and Interactions of Components of the LevQRST Signal Transduction Complex of Streptococcus mutans

Transcription of the genes for a fructan hydrolase (fruA) and a fructose/mannose sugar:phosphotransferase permease (levDEFG) in Streptococcus mutans is activated by a four-component regulatory system consisting of a histidine kinase (LevS), a response regulator (LevR) and two carbohydrate-binding proteins (LevQT). The expression of the fruA and levD operons was at baseline in a levQ mutant and substantially decreased in a levT null mutant, with lower expression with the cognate inducers fructose or mannose, but slightly higher expression in glucose or galactose. A strain expressing levQ with two point mutations (E170A/F292S) did not require inducers to activate gene expression and displayed altered levD expression when growing on various carbohydrates, including cellobiose. Linker-scanning (LS) mutagenesis was used to generate three libraries of mutants of levQ, levS and levT that displayed various levels of altered substrate specificity and of fruA/levD gene expression. The data support that LevQ and LevT are intimately involved in the sensing of carbohydrate signals, and that LevQ appears to be required for the integrity of the signal transduction complex, apparently by interacting with the sensor kinase LevS

Revisiting the pH-gated conformational switch on the activities of HisKA-family histidine kinases

Histidine is a versatile residue playing key roles in enzyme catalysis thanks to the chemistry of its imidazole group that can serve as nucleophile, general acid or base depending on its protonation state. In bacteria, signal transduction relies on two-component systems (TCS) which comprise a sensor histidine kinase (HK) containing a phosphorylatable catalytic His with phosphotransfer and phosphatase activities over an effector response regulator. Recently, a pH-gated model has been postulated to regulate the phosphatase activity of HisKA HKs based on the pH-dependent rotamer switch of the phosphorylatable His. Here, we have revisited this model from a structural and functional perspective on HK853-RR468 and EnvZ-OmpR TCS, the prototypical HisKA HKs. We have found that the rotamer of His is not influenced by the environmental pH, ruling out a pH-gated model and confirming that the chemistry of the His is responsible for the decrease in the phosphatase activity at acidic pH

Haloarchaea swim slowly for optimal chemotactic efficiency in low nutrient environments

Archaea have evolved to survive in some of the most extreme environments on earth. Life in extreme, nutrient-poor conditions gives the opportunity to probe fundamental energy limitations on movement and response to stimuli, two essential markers of living systems. Here we use three-dimensional holographic microscopy and computer simulations to reveal that halophilic archaea achieve chemotaxis with power requirements one hundred-fold lower than common eubacterial model systems. Their swimming direction is stabilised by their flagella (archaella), enhancing directional persistence in a manner similar to that displayed by eubacteria, albeit with a different motility apparatus. Our experiments and simulations reveal that the cells are capable of slow but deterministic chemotaxis up a chemical gradient, in a biased random walk at the thermodynamic limit

The Pseudomonas aeruginosa Chemotaxis Methyltransferase CheR1 Impacts on Bacterial Surface Sampling

The characterization of factors contributing to the formation and development of surface-associated bacterial communities known as biofilms has become an area of intense interest since biofilms have a major impact on human health, the environment and industry. Various studies have demonstrated that motility, including swimming, swarming and twitching, seems to play an important role in the surface colonization and establishment of structured biofilms. Thereby, the impact of chemotaxis on biofilm formation has been less intensively studied. Pseudomonas aeruginosa has a very complex chemosensory system with two Che systems implicated in flagella-mediated motility. In this study, we demonstrate that the chemotaxis protein CheR1 is a methyltransferase that binds S-adenosylmethionine and transfers a methyl group from this methyl donor to the chemoreceptor PctA, an activity which can be stimulated by the attractant serine but not by glutamine. We furthermore demonstrate that CheR1 does not only play a role in flagella-mediated chemotaxis but that its activity is essential for the formation and maintenance of bacterial biofilm structures. We propose a model in which motility and chemotaxis impact on initial attachment processes, dispersion and reattachment and increase the efficiency and frequency of surface sampling in P. aeruginosa